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■ Abstract 

' We apply the no-core shell model to the nuclear structure of odd-mass nuclei straddling 48 Ca. 
Starting with the NN interaction, that fits two-body scattering and bound state data we evaluate 
the nuclear properties of A — 47 and A = 49 nuclei while preserving all the underlying sym- 

■ metries. Due to model space limitations and the absence of 3-body interactions, we incorporate 
phenomenological interaction terms determined by fits to A = 48 nuclei in a previous effort. Our 
modified Hamiltonian produces reasonable spectra for these odd mass nuclei. In addition to the 

i— i differences in single-particle basis states, the absence of a single-particle Hamiltonian in our no- 
t-H ■ core approach complicates comparisons with valence effective NN interactions. We focus on purely 
off-diagonal two-body matrix elements since they are not affected by ambiguities in the different 



(N 



O 



roles for one-body potentials and we compare selected sets of /p-shell matrix elements of our initial 



and modified Hamiltonians in the harmonic oscillator basis with those of a recent model /p-shell 
interaction, the GXPF1 interaction of Honma, Otsuka, Brown and Mizusaki. While some signif- 



icant differences emerge from these comparisons, there is an overall reasonably good correlation 
between our off-diagonal matrix elements and those of GXPF1. 



H I. INTRODUCTION 

_ _ i 

The low-lying levels of the A ~ 47 — 49 nuclei have long been of experimental and theoretical 
interest. On the one hand, extensive experimental information about these nuclei is available 
and, on the other hand, this is a suitable nuclear mass region for developing and testing 
effective /p-shell Hamiltonians. Numerous detailed spectroscopic calculations have been reported. 
For example, in Ref. [3], using a shell model approach, Martinez-Pinedo, Zuker, Poves and Caurier 
have performed full /p-shell calculations for the A = 47 and A = 49 isotopes of Ca, Sc, Ti, 
V, Cr and Mn. They employed the KB3 interaction |8j with phenomenological adjustments and 
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they performed complete diagonalizations to obtain very good agreement with the experimental 
level schemes, transition rates and static moments. Extensive discussions of /p-shell effective 
Hamiltonians and nuclear properties can be found in recent shell model review articles Q-Q]. 

Our interest in these nuclei stems from our goal to extend the no-core shell model (NCSM) 
applications to heavier systems than previously investigated. Until recently, the NCSM, which 
treats all nucleons on an equal footing, had been limited to nuclei up through A — 16. However, 
in a recent paper Q| we reported the first NCSM results for 48 Ca, 48 Sc and 48 Ti isotopes, with 
derived and phenomenological two-body Hamiltonians. These three nuclei are involved in double- 
beta decay of Ca, and the interest in developing nuclear structure models for describing them 
is also related to the need for accurate calculations of the nuclear matrix elements involved in 
this decay. Our first goals were to see the limitations of such an approach applied to heavier 
systems and how much improvement one can obtain by adding phenomenological two-body terms 
involving all nucleons. In brief, the results were the following i) one finds that the charge 
dependence of the bulk binding energy of eight A = 48 nuclei is reasonably described with an 
effective Hamiltonian derived from CD-Bonn inter action [ll in the very limited iV max = basis 
space, while there is an overall underbinding by about 0.4 MeV/nucleon; ii) the resulting spectra 
are too compressed compared with experiment; iii) when isospin-dependent central terms plus a 
tensor interaction are added to the Hamiltonian, one achieves accurate total binding energies for 
eight A — 48 nuclei and reasonable low-lying spectra for the three nuclei involved in double-beta 
decay. Only five input data were used to determine the phenomenological terms - the total binding 
of 48 Ca, 48 Sc, and 48 Ti along with the lowest positive and negative parity excitations of 48 Ca. The 
negative parity calculations are performed in the iV max = 1 basis space. Since the NCSM effective 
2-body interaction is solely responsible for the spectroscopy and involves the interactions of all 48 
nucleons, no single-particle energies are employed. 

In the present paper we extend our previous approach to the odd- A isotopes 47 Ca, 49 Ca, 47 Sc and 
47 K, which differ by one nucleon from 48 Ca. One of our goals is to test whether the same modified 
effective 2-body Hamiltonian used for A — 48 isotopes, is able to describe these odd-A nuclei and 
whether experimentally known single-particle properties emerge in a natural manner. A particular 
feature of the spectroscopy of these odd nuclei is that the spin-orbit splitting gives rise to a sizable 
energy gap in the /p-shell between the f^M and other orbitals (pi/2, Pz/2, fb/2) an d we wanted to 
see if this feature is reproduced in the NCSM. Also, in spite of the differences in frameworks with 
and without a core, we aim to compare selected aspects of our initial and modified Hamiltonian 
with a recent /p-shell interaction, the GXPF1, developed by Honma, Otsuka, Brown and Mizusaki 



121 . We feel it is valuable to compare various /p-shell interactions in order to understand better 



their shortcomings and their regimes of applicability. It is worth mentioning that our interaction 
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and Honma et al. GXPF1 interactions were also tested recently within the framework of spectral 
distribution theory in Ref. [lj| to illustrate their similarities and differences. 

We note that direct comparisons of our matrix elements with those of GXPF1 involve approxi- 
mations that we have attempted to minimize. We specifically list three aspects of the comparisons. 
First, the yl-dependence of GXPF1 is A~ 3 while the A-dependence of the derived effective in- 
teraction we employ is not anticipated to be a simple scale factor. Thus we focus our attention 
on a narrow range of A for the comparison. Second, a more precise comparison would involve the 
derivation of a pure 'valence-only' effective interaction and an appropriate scheme for doing this has 
recently been shown to yield important 3-body forces [ijj . We hope the results presented here 
help motivate that major undertaking. Third, the results we present may also be compared with 
the earlier Brueckner-based matrix elements [17] for the same region since our phenomenological 
adjustments may be expected to simulate the physics accounted for in the perturbative treatment 
of core-polarization. 

Our paper is organized as follows: in Section 2 we give a short review of the NCSM approach and 
we refer the reader to the bibliography for more details. Section 3 is devoted to the presentation of 
our results. Binding energies, excitation spectra, single-particle characteristics, monopole matrix 
elements and matrix element correlations are discussed in subsections along with corresponding 
figures. In the last Section we present the conclusion and the outlook of our work. 

II. NO-CORE SHELL MODEL 



The NCSM 18-23J is based on an effective Hamiltonian derived from realistic "bare" interactions 
and acting within a finite Hilbert space. All A-nucleons are treated on an equal footing. The 
approach is both computationally tractable and demonstrably convergent to the exact result of 
the full (infinite) Hilbert space. 



Initial investigations used two-body interactions 18] based on a G-matrix approach. Later, a 



similarity transformation procedure based on Okubo's pioneering work [24j was implemented to 
derive two-body and three-body effective interactions from realistic NN and NNN interactions. 
Diagonalization and evaluation of observables from effective operators created with the same trans- 
formations are carried out on high-performance parallel computers. 

A. Effective Hamiltonian 



In order to clarify the distinctions from the conventional shell model with a core, we briefly 
outline the ab initio NCSM approach with NN interactions alone and point the reader to the 
literature for the extensions to include NNN interactions. We begin with the purely intrinsic 



4 



Hamiltonian for the A-nucleon system, i.e., 



t<j i<j=i 



where m is the nucleon mass and V^{ri—Tj), the NN interaction, with both strong and electromag- 
netic components. Note the absence of a phenomenological single-particlepotential. We may use 



25| or momentum-space 



either coordinate-space NN potentials, such as the Argonne potentials 
dependent NN potentials, such as the CD-Bonn [ill ]. 

Next, we add to ([1]) the center-of-mass Harmonic Oscillator (HO) Hamiltonian Hqm — 7cm + 
Ucm, where Ucm = \AmVL 2 R 2 , R = \ Y^t=i^i- At convergence, the added -ffcM term has 
no influence on the intrinsic properties. However, when we introduce the cluster approximation 
below, the added -ffcM term facilitates convergence to exact results with increasing basis size. The 



modified Hamiltonian, with pseudo-dependence on the HO frequency fi, can be cast as: 



A r - 1 A r mfi 2 



Kj=l 



(2) 



Next, we introduce a unitary transformation, which is designed to accommodate the short- 
range two-body correlations in a nucleus, by choosing an anti-hermitian operator S, acting only 
on intrinsic coordinates, such that T~L = e^ s H^e s . In this approach, S is determined by the 
requirements that T~L and have the same symmetries and eigenspectra over the subspace K, of the 
full Hilbert space. In general, both S and the transformed Hamiltonian are ^4-body operators. The 
simplest, non-trivial approximation to % is to develop a two-body (a = 2) effective Hamiltonian, 
where the upper bound of the summations "A" is replaced by "a" , but the coefficients remain 
unchanged. 

The full Hilbert space is divided into a finite model space ("P-space") and a complementary 
infinite space ("Q-space"), using the projectors P and Q with P + Q = l. We determine the 
transformation operator S a from the decoupling condition Q a e~ si ' H^e si ] P a — and the simul- 
taneous restrictions P a S^ a 'P a — Q a S^Q a = 0. The a-nucleon-state projectors (P a ,Qa) follow 
from the definitions of the j4-nucleon projectors P, Q. 

In the limit a — ► A, we obtain the exact solutions for dp states of the full problem for any finite 
basis space, with flexibility for choice of physical states subject to certain conditions [2^. This 
approach has a significant residual freedom through an arbitrary residual P a _ space unitary trans- 
formation that leaves the a-cluster properties invariant. Of course, the j4-body results obtained 
with the a-body cluster approximation are not invariant under this residual transformation. It 
may be worthwhile, in a future effort, to exploit this residual freedom to accelerate convergence in 
practical applications. 
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The model space, P2, is defined by AT m via the maximal number of allowed HO quanta of the 



A-nucleon basis states, Nm, where the sum of the nucleons' 2n + I < N n 



N 



spsmin 



N M , and 



where N m 



denotes the minimal possible HO quanta of the spectators, nucleons not involved in 



the interaction. For example, in °B we have N spsm i n = 4 as there are 6 nucleons in the Op-shell 
in the lowest HO configuration and N m = 2 + iV max , where iV max represents the maximum HO 
quanta of the many-body excitation above the unperturbed ground-state configuration. For 10 B 
with A^m = 12, N m — 8 for an iV max = 6 or %hft" calculation. With the cluster approximation, a 
dependence of our results on N max (or equivalently, on N m or Nyi) and on Q arises. The residual 
JV max and fl dependences will infer the uncertainty in our results arising from effects associated 
with increasing a and/or effects with increasing iV max . In the present work, we retain the iV max = 
basis space and hfl = lOMeV employed in Ref. [9(. 

At this stage we also add the term Hcm again with a large positive coefficient (the Lagrange 
multiplier) to separate the physically interesting states with Os CM motion from those with excited 
CM motion. We diagonalize the effective Hamiltonian with the m-scheme Lanczos method to obtain 
the P-space eigenvalues and eigenvectors 10] . All observables are then evaluated free of CM motion 



effects 



10j. In principle, all observables require the same transformation as implemented for the 



Hamiltonian. We obtain small renormalization effects on long range operators such as the rms 
radius operator and the B(E2) operator when we transform them to P-space effective operators 



at the a = 2 cluster level 



22, 



28|. On the other hand, when a = 2, substantial renormalization was 
observed for the kinetic energy operator [27]. and for higher momentum transfer observables j^ . 
Recent applications and extensions include: 



(a) spectra and transition rates in p-shell nuclei 

(b) comparisons between NCSM and Hartree-Fock [30I ] : 

(c) neutrino cross sections on 12 C |31|; 

(d) novel AW interactions using inverse scattering theory plus NCSM 

(e) solving for light nuclei with chiral effective interactions 

(f ) studies of alternative converging sequences [36( ; 
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34|; 
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37]; 



(g) solving for light nuclei with alternative renormalization schemes 

(h) development of effective interactions [15] and operators [3] with a core; 

(i) solving for light nuclei with bare AW interactions with extrapolation to the infinite matrix 
limit 381 ; 



G 



(j) phenomenology of hadronic structure and quantum field theory 



39] 



We close this theory overview by referring to the added phenomenological NN interaction terms 
found adequate for obtaining good descriptions of A = 48 nuclei jsj. Three terms are added - 
central modified Gaussians with isospin dependent strengths and a tensor term. The forms were 
chosen primarily for their simplicity. However, based on experience in light nuclei, it is desirable to 
include a tensor term in order to better fit the spin-sensitive properties. In the hope of obtaining 
a simple NCSM Hamiltonian for the BE and spectra of A = 47, 48 and 49 nuclei, we assume the 
added NN interaction terms are A-independent. That is, they are taken to complement the A- 
dependent NCSM H e ff in a very simple manner. NCSM results obtained below with the modified 
Hamiltonian are referred to with "CD-Bonn + 3 terms" results. It is our hope that these terms 
accommodate, to a large extent, the missing many-body forces, both real and effective. This hope 
will be tested in the future when increasing computational resources will allow larger basis spaces, 
improved a — 3 and a = 4 calculations as well as the introduction of true NNN and NNNN 
potentials. 



III. RESULTS AND DISCUSSION 
A. Binding energies 



First we present the calculated total interaction energies (Hamiltonian ground state eigenvalues) 
in Fig. Q] which we compare with experiment. One observes that ground states calculated with 
our derived ab initio H e f f lie above the experimental values by approximately 20 MeV. This shift 
is similar to that observed in the case of all A = 48 isotopes [9j. We note that with CD-Bonn we 
have nearly the same increase in binding from 47 Ca to 48 Ca as from 48 Ca to 49 Ca, which signals a 
lack of subshell closure. 

For the modified Hamiltonian (CD-Bonn + 3 terms) the NCSM produces reasonable agreement 
with experiment with deviations much less than 1% as seen in Fig. [TJ There is a simple spreading 
of the theoretical ground states relative to experiment. In particular, we now observe the desired 
subshell closure condition where the increased binding from 47 Ca to 48 Ca significantly exceeds that 
from 48 Ca to 49 Ca. 



B. Excitation energy spectra 

The excitation energy spectra for 49 Ca, 47 Ca, 49 Sc and 47 K are shown in Figs. OS] respectively. In 
every case the ab initio NCSM results with CD-Bonn are far too compressed relative to experiment 
- a feature also seen in the A = 48 results Here, we trace this primary defect to the inferred 
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FIG. 1: The experimental and theoretical ground state 
energy levels for A — 47 — 49. The results in the second 
and third columns are labelled by their Hamiltonians. 
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FIG. 2: Experimental and theoretical excitation energy lev- 
els for 49 Co. Both CD-Bonn and CD-Bonn+3terms results 
are presented. 



properties of the neutron orbits. That is, the incorrect ground state spin seen in Fig. [5] and the 
absence of a significant excitation energy gap in Fig. [3] indicate the spin-orbit splitting of the 
neutrons is insufficient to provide proper subshell closure at the neutron O/7/2 orbit. This defect 
is rectified in the results with CD-Bonn + 3 terms as seen by the corresponding spectra in Figs. [2] 
and [3J Similar tendencies have been seen before with valence G- matrix interactions and identified 
as a problem with the L 2 dependence of the single-particle states 

The CD-Bonn results in Figs. |4]and[5]are more difficult to interpret due to the glaring deficiencies 
just mentioned for the neutrons with the CD-Bonn Hamiltonian. We will show below that the 
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FIG. 3: Experimental and theoretical excitation en- 
ergy levels for 47 Ca. Both CD-Bonn and CD- 
Bonn+3terms results are presented. 




FIG. 4: Experimental and theoretical excitation 
energy levels for 49 Sc. Both CD-Bonn and CD- 
Bonn+3terms results are presented. 



proton shell closure is better established with CD-Bonn. This supports the assertion that the 
main deficiencies seen in the third columns of Figs. 2] and [5] are indeed likely to reside with the 
inferred neutron spin-orbit splitting problem. 

The modified Hamiltonian provides greatly improved spectra for all four nuclei as seen in the 
second columns of Figs. [2][5j It is to be noted that these nuclei were not involved in the fitting 
procedure used to determine the parameters of the added phenomenological terms. Perhaps the 
most significant remaining deficiency is the incorrect ground state spin for K as seen in Fig. [5] 
This is the first case of a nucleus in the region of A — 47 to A — 49 (12 nuclei studied to date) 
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FIG. 5: Experimental and theoretical excitation 
energy levels for 47 K. Both CD-Bonn and CD- 
Bonn+3terms results are presented. 
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FIG. 6: Experimental and theoretical levels that are 
dominantly single-proton and single-neutron particles 
or holes of 4S Ca. The levels are labeled by (n, 1, 2j), 
and the dashed lines are the Fermi energies. 



where we did not obtain the correct ground state spin with CD-Bonn + 3 terms Hamiltonian. 



C. Single-particle characteristics 



In order to better understand the underlying physics of our NCSM results, we investigate the 
single-particle-like properties of our solutions. 

In a simple closed-shell nucleus, we expect the leading configuration of the ground state solution 
in our m-scheme treatment to be a single Slater determinant. Single-particle (or hole) excitations 
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should be easily identified by the character of their leading configurations - i.e. a single-particle 
creation (or destruction) operator acting on the ground state Slater determinant of the reference 
nucleus 48 Ca. For our odd-mass nuclei, this is the character we seek. That is, we take the standard 
phenomenological shell model configuration of a single Slater determinant with a closed sd-shell 
for the protons and a closed f^/ 2 subshell for the neutrons and look for the appropriate states 
which have a single nucleon added to (or subtracted from) that Slater determinant. We accept 
states as "single-particle-like" when we find one with a leading configuration having more than 
50% probability to be in the simple configuration just described. When the majority weight is 
distributed over a few states, we use the centroid and we discuss those cases in some detail below. 
We were not successful in locating all the expected single-particle-like and single-hole-like states. 
That is, those absent from our presentation below were spread among a large number of eigenstates. 

For a closed-shell nucleus (Z, N) the single-particle energies (SPE) for states above the Fermi 
surface are related to the binding energy differences: 

e> =BE(Z,N) -BE*(Z+1,N) 
and 

e> = BE(Z,N) - BE*(Z,N + 1). 

The SPE for sates below the Fermi surface are given by 
e< = BE*(Z -1,N) - BE(Z,N) 
and 

e< = BE*(Z,N- 1) -BE(Z,N). 

The BE are ground state binding energies which are taken as positive values, and e will be 
negative for bound states. (BE* = BE — E x ) is the ground state binding energy minus the 
excitation energy of the excited states associated with the single-particle states. 

Experimental SPE's and the results of our analysis are shown in Fig. [51 The experimental SPE's 
for protons and neutrons follow B.A.Brown's analysis [4|. To guide the eye, we draw a horizontal 
line to indicate the vicinity of the Fermi surfaces for the protons and neutrons. 

Fig. [6] shows that proton shell closure is established with both Hamiltonians, the CD-Bonn and 
the CD-Bonn+3 terms. The correct energy locations are better approximated with the modified 
Hamiltonian. Fig. [6] also shows that neutron subshell closure only appears with the modified 
Hamiltonian. Here, the ordering is correct but the states are considerably more spread out com- 
pared with experiment. 

Let us consider some of the details underlying the single-particle-like states. The situation for 
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the IP3/2 or "lp3" state in the left panel of Fig. |6l the proton single-particle state in 49 Sc with 
the modified Hamiltonian, is quite interesting. It appears that this state is mixed over several 
excited states in the spectrum. We can take the strength spread over several states and construct 
a centroid for this lp3 state by a weighted average over the states carrying that strength. Here are 
the relevant input ingredients. 

The first excited state of 49 Sc is a 3/2 , as seen in the second column of Fig. HJ with about 51% 
of the occupancy of the lp3 state. Its eigenvalue is —425.151 MeV compared to a ground state of 
—428.365 MeV. The 18th state in the 49 5c spectrum is also a 3/2~ with 28% of the occupancy 
of the lp3 state. Its eigenvalue is -422.803 MeV. The 24th state is also a 3/2~ with 21% of the 
occupancy of the lp3 state. Its eigenvalue is —422.440 MeV. 

Thus, to a good approximation, the lp3 strength is spread over these three states. We will 
identify the weighted average [0.51 x (-425.151) +0.28 x (-422.803) +0.21 x (-422.440)] = -423.79 
as the centroid of the single particle lp3 state which we then include accordingly in the second 
column of the figure. 

For the proton hole states with the modified Hamiltonian, we perform a detailed search up to 
excitation energies of about 14 MeV in the 47 K spectra. It appears that the 0d 5 / 2 single- hole state 
is spread among many states with the largest observed concentration on the 5/2 + state at —386.17 
MeV (13.36 MeV of excitation energy). Here, we find a single J v = 5/2 + state in 47 K with 30 % 
0^5/2 vacancy and we assign this state to our Qd 5 / 2 single-hole state. Most of the 0d 5 / 2 strength, 
however, was not observed among the limited number of converged eigenstates. 

Let us consider the 49 Ca results with the modified Hamiltonian in the upper right panel of Fig. 
O The ground state is approximately a pure [(IP3/2) 1 (0f7/2) 8 ) configuration. We note that the 
spacing for the subshell closure is in good agreement with experiment while there is a shift of a 
couple MeV towards more binding in the model as previously indicated in Fig. [TJ A nearly pure 
lpi/2 single-particle state is obtained at 5.235 MeV excitation energy and an extra low- lying 7/2~ 
appears with 2p — lh character (see Fig. [2]). Our lowest-lying 5/2~ consists of 2p — lh character 
relative to subshell closure. 

We contrast the modified Hamiltonian's results for the 49 Ca ground state with those obtained 
using the ab initio CD-Bonn where [(lp3/2) 4 (lpi/2) 2 (0/7/2) 3 ]i/2- is the dominant configuration 
reflecting again the inadequacies of the neutron single-particle properties with CD-Bonn. 



12 



D. Monopole matrix elements V(ab;T) 



The monopole matrix element is defined by an angular momentum average of coupled doubly- 
reduced two-body matrix elements: 

vt h T\ HA2J+l)V(abab; JT) 

V(ab;T) = S J (2J+1) ' (3) 

For our NCSM Hamiltonians the "V appearing in Eqn. 

[3] signifies the full 2-body intrinsic-coordinate Hamiltonian, T re i + V e ff, except that we omit the 
Coulomb interaction from this analysis. 

We examined the monopole character of our initial CD-Bonn Hamiltonian and we compared it 
with the monopole features of the GXPF1 interaction [12j. Although the monopole characteristics 
were similar, we do not present the detailed comparisons here due to the ambiguity of the role 
of the SPE's. That is, one may shift some Hamiltonian components between SPE's and two- 
body matrix elements (TBME's) and this obscures direct comparisons of a subset of our TBME's 
with the corresponding subset of GXPF1. In order to summarize a comparison of the underlying 
theoretical interactions, we list in Table 1 a simplified overview of their differences and similarities. 

For a sample comparison of the interactions, we present a small set of two-body /p-shell matrix 
elements applicable to the present investigation in Table 2. For convenience we present two columns 
of key differences in the matrix elements: "diffl" represents the difference between the G-matrix 
and the GXPF1 interaction resulting from adjusting the G-matrix elements to fit spectra; and, 
"diff2" represents the difference between our ab initio H e ff and our modified H e ff. The scale of 
the changes from the respective starting solutions appears comparable, though one of the "diff2" 
values reaches -1.2123 MeV. 



E. Matrix element correlations 



We present in Figs. [JJ -QJ] the correlations between pairs of /p-shell interaction matrix element 
sets. With Fig. [JJ we observe the high degree of correlation between the 195 matrix elements of 
our starting Hamiltonian, CD-Bonn, and our modified Hamiltonian, CD-Bonn + 3 terms. This 
indicates that, for the most part, our Hamiltonian is minimally modified by the addition of the 
phenomenological terms. Such a high correlation is reminiscent of the high correlations seen 
between GXPF1 and its starting interaction, the G-matrix [ijj]. It is interesting to see if certain 
groups of matrix elements appear to be more correlated than the others. We distinguish the 
diagonal V(abab; JT) matrix elements that contribute to the monopole by different symbols. Filled 
circles stand for V(abab; JT = 0) and filled squares stand for V(abab\ JT = 1) matrix elements. 
All the remaining matrix elements V(abcd; JT) where at least one single-particle-state (sps) of 



13 



the bra is different from a sps of the ket are plotted as open circles for T = and open squares 
for T = 1. We see the filled square points, that correspond to V(abab; JT =1) matrix elements, 
are farther from the linear fit, ranging between 1 MeV and 2 MeV away from the linear fit line. 
Therefore, these monopole matrix elements have received larger corrections than others in the 
process of fitting the A = 48 isotopes. 

In the forth paragraph of the same subsection we changed the sentence "The colors and the 
symbols used are the same as those in figures 7 and 8." with "The symbols used are the same as 
those in figures 7 and 8." 

To see the stronger correlations more clearly, we next choose in Fig. [8] to eliminate all matrix 
elements contributing to the monopole. There are 135 remaining matrix elements out of the 195 
total. The degree of correlation between the 135 matrix elements significantly improves with much 
less deviation from the linear fit. We can see another feature of the correlation by comparing the 
linear fit with the 45-degree line in Fig. [8] (similar pattern seen in Fig. [7]). On the one hand, 
we see that the CD-Bonn + 3 terms matrix elements are shifted towards greater attraction where 
CD-Bonn is already attractive. On the other hand, the CD-Bonn + 3 terms matrix elements are 
shifted towards greater repulsion where CD-Bonn is already repulsive. Overall, we observe that the 
larger differences between the CD-Bonn and CD-Bonn + 3 terms are coming from the monopole 
terms. This seems natural in light of the fact that the phenomenological terms have the effect 
of adjusting the single particle features of the theory towards agreement with experiment. That 
is, the monopole terms receive the largest adjustments as required to achieve the needed single 
particle features. 

It is then very interesting to observe in Fig. |H] the lack of correlation between our starting 
Hamiltonian, CD-Bonn, and the G-matrix underlying the GXPF1 interaction. Points are generally 
farther away from the fit line than in the correlation of CD-Bonn + 3 terms with CD-Bonn case. 
Note that the G-matrix is a renormalization procedure and the specific results for GXPF1 are 
developed from the bare CD-Bonn interaction. This lack of correlation in Fig. [9] reflects the major 
differences in the underlying effective interaction theories that are summarized in Table 1. 

We now make the same set of comparisons between CD-Bonn and G-matrix in Figs. [9] and 
ITTJ1 as we performed in Figs. [7J and [5] The symbols used are the same as those in Figs. [7J and 
[8] Comparing Fig. [9] and Fig. [10] we can see how, after eliminating the matrix elements that 
contribute to the monopole, the correlation is significantly improved with the linear fit in Fig. [10] 
and it now overlaps well with the 45-degree line. Again, this shows similarities between these two 
interactions with the difference arising primarily from the monopole part. 

We can comment further about the comparison presented in Fig. [9] by observing that the full 
Hamiltonian developed from the G-matrix includes single-particle energy (SPE) contributions. On 
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the other hand CD-Bonn does not have additional SPE contributions since those contributions 
are already included in the 2-body matrix elements. In fact, those SPE contributions to our 
interactions are embedded in the monopole terms. This is one reason why the correlations improve 
when we proceed from Fig. [S]to Fig. I1Q[ removing the monopole terms. 

Finally, in order to focus as clearly as possible on the 2-body interaction effects , we present in 
Fig. [IT] the correlation of matrix elements V(abcd; JT)(A = 48) between CD-Bonn+3terms and 
GXPF1, where we retain only those that cannot contribute to a single-particle Hamiltonian. That 
is, we eliminate all two-body matrix elements where at least one single-particle-state (sps) of the 
bra equals a sps of the ket. There are 56 remaining two-body matrix elements. Differences ranging 
up to about 3 MeV are observed which should lead to differences in experimental observables. 
Comparisons of spectra and other properties with these Hamiltonians, as one proceeds further 
from A = 48, could shed more light on their differences. 

IV. CONCLUSIONS AND OUTLOOK 

We have presented an initial NCSM investigation of the spectral properties of the A = 47 and 
A = 49 nuclei that are one nucleon away from doubly-magic 48 Ca. We have shown that the NCSM 
with a previously introduced modified Hamiltonian produces spectral properties in reasonable 
accord with experiment. Shell closure and single-particle spectral properties are obtained indicating 
a path has been opened for multi-shell investigations of these nuclei within the NCSM. We are 
undertaking such additional investigations. Also, for a better understanding of various /p-shell 
interactions we made a comparison between our initial and modified /p-shell matrix elements in the 
harmonic oscillator basis with the GXPF1 interaction [12J. Our initial and modified NCSM H e ff 
matrix elements in the /p-shell are strongly correlated. We found some evidence suggesting that 
significant differences in single-particle properties may underly some of the distinctions between 
our H e ff and the GXPF1 interaction. The differences were reduced when we compared the purely 
off-diagonal matrix elements in Fig. 111! Additional applications could reveal the importance of 
these distinctions in greater detail. 
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Hamiltonian Property 


G-matrix 


NCSM cluster H eff 


Oscillator parameter dependence 


Yes 


Yes 


Depends on the choice of P-space 


Yes 


Yes 


Reguires effective multi-nucleon 






interactions as corrections 


Yes 


Yes 


Translationally invariant 


No 


Yes 


Starting energy dependence 


Yes 


No 


Single-particle spectra dependence 


Yes 


No 


A-dependence 


No 


Yes 



TABLE I: Overview of the differences and similarities of the two theoretical approaches that underlie the 
Hamiltonians whose matrix elements are compared in this work. 



2ja 2j b 


2jc 2j d 


J T 


G 


GXPF1 


diffl 


CD-Bonn CD-Bonn 
+3 terms 


diff2 


7 


3 


7 


3 


5 





-2.1167 


-2.8504 


-0.7337 


-1.0390 


-1.3413 


-0.3023 


3 


3 


5 


5 





1 


-0.5243 


-1.1968 


-0.6725 


-0.6019 


-1.1129 


-0.5109 


7 


7 


7 


7 


3 





-0.2309 


-0.8087 


-0.5778 


0.5597 


0.5555 


-0.0042 


7 


5 


7 


5 


6 





-2.3465 


-2.9159 


-0.5693 


-1.3743 


-1.8599 


-0.4856 


7 


5 


7 


5 


5 





-0.0203 


-0.5845 


-0.5642 


0.5813 


0.4117 


-0.1693 


3 


1 


3 


1 


2 


1 


-0.7965 


-0.2822 


0.5143 


-0.0068 


-0.4932 


-0.4864 


7 


7 


5 


5 





1 


-1.9095 


-1.3288 


0.5806 


-2.2586 


-3.4709 


-1.2123 



TABLE II: Comparison of selected two-body matrix elements V(abcd; JT) (Mev) (A = 48) for which the 
difference between our interaction is large, "diffl" represents the difference between GXPF1 and G while 
"diff2" is the difference between CD-Bonn+3terms and CD-Bonn. 
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1 (A=48) 



FIG. 7: (color online) Correlation of V(abcd; JT) matrix 
elements between CD-Bonn+3terms and CD-Bonn. The 
matrix elements of T = and T = 1 are shown by circles 
and squares, respectively. The filled circles and the filled 
squares are for all V(abab; JT) matrix elements that con- 
tribute to the monopole V(ab;T). The open circles and 
open squares are for all the remaining matrix elements. 
There are no monopole shifts. The solid straight line 
represents a linear fit to all the matrix elements. The 
diagonal dashed line represents the reference correlation 
line at 45-degrees. 
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FIG. 8: (color online) Correlation of V(abcd;JT) 
between CD-Bonn+3terms and CD-Bonn, where we 
removed all V(abab; JT) matrix elements that con- 
tribute to the monopole V(ab;T). The solid straight 
line represents a linear fit to all the plotted points. 
The open circles and the open squares stand for T — 
and T — 1, respectively. The diagonal dashed line 
represents the reference correlation line at 45-degrees. 




FIG. 10: (color online) Correlation of V(abcd; JT) 
between CD-Bonn and G. See the caption to Fig. [8] 



19 



3 
2 



1 1 1 1 


I L 

/ / ■ 




/ / 


- 






° / / 




x " 

F ■ 






- / X n n 








/ / 








s / " 




7 / / 




/ / 




■_ / / 




: / / 




./ / 

■ i i i 





-5 -4-3-2-10123 

GXPF1 



FIG. 11: (color online) Correlation of the matrix el- 
ements V(abcd; JT) between CD-Bonn+3terms and 
GXPF1, where we retain only the off-diagonal ma- 
trix elements, i.e. the 56 matrix elements that can- 
not contribute to a single-particle Hamiltonian (see 
text). The open circles and the open squares stand 
for T = and T = 1, respectively. The thick solid line 
represents a linear fit to all the plotted points. The 
diagonal dashed line represents the reference correla- 
tion line at 45-degrees. 



